Josephson junction qubit network with current-controlled interaction 
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We design and evaluate a scalable charge qubit chain network with controllable current-current 
coupling of neighbouring qubit loops via local dc-current gates. The network allows construction of 
general N-qubit gates. The proposed design is in line with current main stream experiments 



Although a working solid-state quantum computer 
with hundreds of qubits remains a distant goal, coupling 
of a few solid-state qubits is now becoming feasible. Sev- 
eral groups have succeeded to operate single Josephson 
junction (J J) qubits [1-8], but the art of multiple J J qubit 
gates is still in its infancy. A few challenging experiments 
with coupled JJ-qubits have been reported [9-13]. How- 
ever, so far experiments on coupled JJ qubits have been 
performed without direct physical control of the qubit- 
qubit coupling. 

There are many proposed schemes for two(multi)-qubit 
gates where an effective qubit coupling is controlled by 
tunings of qubits or bus resonators [14-17]. However, 
there are also suggestions how to control physical qubit 
interaction [19-23], most of which require local magnetic 
field control. Recently, Yamamoto et at [12] successfully 
implemented a CNOT gate using fixed capacitive cou- 
pling between two charge qubits, controlling the effective 
qubit-qubit interaction by tuning single-qubit level split- 
tings into resonance - however, this method might not 
be well suited for more advanced gates on charge qubits 
because of strong decoherence when qubits are operated 
away from the degeneracy points. 

In this paper we present a new solution for control- 
lable physical qubit-qubit coupling, as shown in Fig. 1. 
The network has the following properties: (a) nearest- 
neighbour qubit-qubit coupling controlled by external 
bias current, (b) qubits parked at the degeneracy points, 
also during qubit-qubit interaction, (c) separate knobs for 
controlling individual qubits and qubit-qubit coupling, 
(d) scalability. An important feature is that the net- 
work is easily fabricated, and is in line with current main- 
stream experiments. 

The network under consideration consists of a chain of 
charge qubits - Single Cooper Pair Transistors (SCPT) 
- with loop-shaped electrodes coupled together by cur- 
rent biased coupling JJs at the loop intersections (Fig. 
1). The loop-shaped electrode was introduced [15,1] to 
provide external control of the Josephson coupling of the 
qubit island to the reservoir. The loop design creates 
an (inductive) interface to the qubit by means of cir- 
culating currents [24] , which has been used as a tool for 
qubit readout by Vion et al. [3]. We employ these current 
states in the qubit loops to create controllable coupling 
of neighbouring qubits. The results of this paper are de- 



rived in the charge qubit limit Ec 3> Ej. However, the 
analysis and the coupling mechanism also apply to the 
case of Ec ~ Ej, describing the charge-phase qubit [3]. 




FIG. 1. Network of loop-shaped SCPT charge qubits, cou- 
pled by large Josephson junctions. The interaction of the 
qubits (i) and (i+ 1) is controlled by the current bias l\ or by 
simultaneous current biasing of readout junctions. Individual 
qubits are controlled by voltage gates, V gi . Single-qubit read- 
out is performed by applying an ac [25, 26] or dc [3] current 
pulse JJ™ to a particular JJ readout junction. Alternatively, 
readout of charge may be performed, eg. using a RF-SET ca- 
pacitively coupled to the island [27] The two Josephson junc- 
tions of the i-th SCPT are assumed to have identical Joseph- 
son energies Eji. The Josephson energies of the coupling JJs 
Ej and the readout JJs Ej 1 (identical for simplicity) are much 
larger than the corresponding charging energies, Ej S> E c . 
<j>\ (4>T) i s the phase difference across the i-th coupling (read- 



out) JJ. For an open N-qubit chain we choose (f>o ~ 4"% 



while for a closed chain 



bo = • 



= 0, 



Left without any external current biasing of the cou- 
pling and readout JJs, the network acts as a quan- 
tum memory of independent qubits (neglecting a weak 
residual interaction, to be discussed below). When a 
bias current is sent through the coupling J J in Fig.l, 
the current-current interaction between the neighbour- 
ing qubits is switched on and increases with increasing 
bias current. Moreover, if both of the readout JJs of the 
same qubits are biased well below threshold, again there 
is nearest-neighbour coupling via the circulating currents 
[29]. However, if one of the readout JJs is current biased, 
this only affects that particular qubit and allows read- 
out of individual qubits. Thus the bias currents serve as 
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the interaction control knobs. The loop inductances are 
assumed to be sufficiently small to neglect qubit-qubit 
coupling via induced magnetic flux, as well as undesir- 
able qubit coupling to the magnetic environment. In ad- 
dition, we assume negligible capacitive coupling between 
the islands, which are well shielded by the injection leads. 

The SCPT qubit chain system shown in Fig. 1 is de- 
scribed by the Hamiltonian 

i 

where H scptA = (E C i/2)(l - n gi )a zi ~ E J{ cosQ l a xl , and 
where H b ?sc - = Q\ /2C b - E^costf - f-J b 4>\ is the 
Hamiltonian of the coupling J J. H™ sc i is the similar 
Hamiltonian of the readout junction, which for simplic- 
ity is chosen with the same parameters. en gi = C g iV g i 
is the induced gate charge on the i-th SCPT island, the 
charging energy of which is defined by Ea = (2e) 2 /2Cs. 
Finally, Qi is the charge on each coupling JJ obeying 
the commutation relations [Qi,<f> b ] = 2ieSij. H SCPTt i has 
been truncated to the two lowest charge states, assuming 
Ea 3> Eji and n gi w 1, and a small correction to the 
charging energy of the coupling JJs (~ Ci/C b <C 1) has 
been neglected. The flux quantization tfin+(f)i2 — 29i = 0, 
in every loop (assuming zero external flux) introduces 
a dependence of the qubit Joscphson energy on the 
phase differences across the coupling and readout JJs, 
9i = (4>\ - 4>\_ L - 4>T)I 2 - Thi s qubit-oscillator interac- 
tion is the origin of the qubit-qubit interaction. 

For proper functioning of the network, the critical con- 
ditions are E b » fkv p ,Ej, where ui p is the plasma fre- 
quency of the coupling JJs, establishing the phase regime 
for the coupling JJs with small fluctuations of phase 
Si = <fi b — 4> b , < Sf > ~ hujp/Ej <C 1, around energy min- 
ima determined by the control current, sm(f) b = I b /I b . 
We only consider the regime of negligible macroscopic 
quantum tunneling (MQT). 

Using a harmonic approximation for the periodic po- 
tential terms in Eq. (1), all coupling JJs are reduced to 

harmonic oscillators with level spacing Uui p = ^j2e b E b J , 

where e b = E b j cos 4> b . Each SCPT term in Eq. (1) is 
then, in the lowest approximation with respect to har- 
monic amplitudes, reduced to a qubit Hamiltonian, 

Hqi = -^i 1 - n g i)a zi - EjiCos9ia x i, (2) 

plus a linear oscillator-qubit interaction, = Xi(Si — 

Si-i)a x i, proportional to the coupling strength A» = 
Eji svnOi and to the phase deviation in the coupling JJs. 
This generates controllable nearest-neighbour qubit in- 
teraction terms which appear only in the presence of 
bias currents and describe displacement of the oscilla- 
tors driven by the qubits. There are also quadratic terms, 
H il]i = (e J /2)(^-(5 J _i) 2 (T ;l;i , where e l = E M cos^, which 



induce relatively small permanent residual qubit coupling 
due to oscillator squeezing driven by the qubits. 

The harmonic oscillators can with good accuracy be as- 
sumed to stay in the ground state during all quantum op- 
erations on the network at low temperature. For current 
control pulse durations T 3> w~ 1 , the probability to ex- 
cite the oscillator due to qubit flips away from the degen- 
eracy point is estimated as W ~ Eq^I — n g i) 2 /Tuxj p e b j <C 
1, when the linear qubit-oscillator coupling is switched on 
(Xi ~ Eji). In the residual interaction regime (Aj = 0), 
the excitation probability is several orders of magnitude 
(Uuip/Ej <C 1) smaller. Hence, we average over the 
ground state of the oscillators and finally arrive at the 
effective Hamiltonian for the qubit network, 

~H = *^2(H ql + rji <7 X i<y x{i+i)) + ^ Kij <Jxi<Jx 3 (3) 

i i=jLj 

where rji = AjAj+i/e* and the energies of the 

controllable and the residual interactions, respectively. 
The maximum controllable interaction energy is a factor 
of Ej/e b smaller than the qubit level splitting, <~ 2Ej. 
The residual qubit-qubit interaction effectively connects 
all of the qubits but it is smaller than the controllable 
coupling by a factor of hu) p /E b -C 1. 

The interaction energy r\i can be expressed in terms of 
the currents Ii — (e/K)Eji smOi circulating in neighbour- 
ing qubit loops as rji = L aft IiI i+ i, where L e{{ = Ti 2 / (4e 2 e^) 
is the effective inductance introduced by the coupling J J. 

In order to exclusively couple the qubits (i) and (i + 1) 
one should apply a nonzero current bias I b , while I b ±1 = 
and I b ± 2 =0. In this case the coupling amplitude is 
given by the equation, 

E JiEj{i+}) 
' AE b j cos $ 2 1 ; 

The coupling is quadratic, ~ (I b /I b ) 2 , for small cur- 
rent bias, and diverges when approaching the critical 
current. An alternative way to switch on the qubit- 
qubit coupling is to apply dc bias currents (below the 
critical value) simultaneously to both of the two neigh- 
boring readout junctions, instead of activating the cou- 
pling junction, resulting in an interaction energy rji = 
-{E M E, ]{l+1) /±E b ) sin(^/2) sin(0j_i/2). 

The present strategy is to park the qubits at the de- 
generacy points, where the coherence time is maximum 
[3], and then to operate with (a) short dc- voltage pulses 
or, alternatively, microwave resonant excitation, to per- 
form singlc-qubit operations with qubit-qubit coupling 
switched off (rj = 0), and with (b) dc-currcnt pulses 
(77 7^ 0) to perform two-qubit rotations at the degeneracy 
points. 

The readout of individual qubits can be performed 
by probing the corresponding junction with ac current 
[25,26], while keeping zero bias at the coupling junctions. 
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Since amplitude of the phase oscillation in the readout 
junction remains small during measurement, our theory 
applies, and neighboring qubits will not be disturbed. 
Another option would be to pulse the dc current through 
the readout junction above the critical value [3]. Re- 
quired that the qubits are operated in the charge regime 
(Ec 3> Ej), readout of the islands charge state is also 
possible by means of a capacitive probe, eg. using a SET 
electrometer. 

We now focus on two-qubit gate operation. The Hamil- 
tonian in Eq. (3) is diagonal in the current basis when 
all qubits are put at the degeneracy points. Considering 
two neighbouring qubits, 1 and 2 in the current eigenbasis 
1 00} ,|01),|10),|11), the Hamiltonian explicitly becomes 
(assuming for simplicity identical qubits e = e\ = e 2 ), 



/2e 



H = 



V 



o 

-m 
o 
o 



o 
o 

-m 
o 



\ 



(5) 



-2e + m J 



where rji is given by Eq. (4), We can now use the cur- 
rent control bias to perform coupled-qubit phase rota- 
tions [16]. We define a basic entangling two-qubit gate 
operation, the — tt/2 zz-rotation, 



. J t dtn(t) 
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(6) 



by choosing the appropriate amplitude I (i.e. r/i) and 
duration T of the bias current pulse, determined by the 
simple integral equations f T (2e/fi)dt — (mod tt) and 
J T (r]i/h) dt = — 7r/4. The operation is only slightly more 
complicated for non-identical qubits. The current pulse 
shape is of no importance, except that it must be adia- 
batic with respect to the harmonic degrees of freedom, 
rji <C hujp. 

By means of the — tt/2 zz-rotation and single qubit 
rotations it is now straightforward to construct any de- 
sired quantum operations, including generalized quan- 
tum gates. Standard two-qubit gates such as the CNOT 
operation require a short sequence of additional single 
qubit tt/2 rotations 



hM z HhP 

>] 



where [k]i = exp[— i(ir/4)Tki], Tki are Pauli matrices in 
the current basis, and the Hadamard operation [H] cor- 
responds to the sequence [x] [z] [x] . Another useful opera- 
tion, CNOT-SWAP, can be also introduced, 



xHz> 



which allows effective implementation of quantum algo- 
rithms on qubit networks with nearest neighbour interac- 
tion [31,33]. The operations have been optimized in the 
sense that the [z] rotations can be performed using the 
natural precession of the qubits. 




FIG. 2. Maximum coupling energy r\ of two neighbouring 
qubits with a current bias applied to the coupling JJ which 
corresponds to the MQT rate Tmqt = (2Ej/h) x 1(T 4 , for 
different shunt capacitances (solid), C b hunt = (a), IOC 6 (6), 
40C 6 (c). The corresponding plasma frequency, for the same 
parameters, (dashed) decreases with increased C h shunt . Inset: 
the ratio of the controllable interaction and the residual in- 
teraction. 

The time needed for a two qubit operation is given by 
the coupling strength 77, whose upper limit is set by MQT 
in coupling JJs, and depends on the plasma frequency: 
a lower plasma frequency yields lower rate of MQT and 
thus higher maximum current bias. Thus, a stronger 
controllable coupling is achieved by adding a large shunt 
capacitance C b shunt to reduce the plasma frequency. On 
the other hand, the latter is restricted by the adiabaticity 
condition, rji <C Tilj p . The maximum coupling energy can 
be estimated using the standard expression describing 
MQT in current biased JJs [30] , neglecting the small cir- 
culating currents, r MQTi j w uip^Os/n) 1 / 2 exp(— s), where 
s = (24^/5to p )(cos^) 5 / 2 and u p = {2E h J E h c ) 1 / 2 is 
the bare plasma frequency of the capacitively shunted 
coupling junction, E b c = (2e) 2 /2{C b + C b shunt ). Re- 
quiring the lower bound for cos6*i to be larger than 
(u>p/E b ) 2 / 5 under condition that the MQT rate remains 
negligibly small, the adiabaticity condition gives, u> p 3> 
Ej(Ej j Ej) 1 ' 4 , and the maximum coupling energy in Eq. 
(4) becomes r\ ma ,x ~ Ej(Ej / Ej) 1 / 2 . Quantitative results 
for the maximum coupling energy are shown in Fig. 2. 
Note, that the residual interaction is reduced for small 
plasma frequency (see inset in Fig. 2). Taking the inter- 



3 



action energy to be rj max = (2_Ej)/100, the time needed 
for the —it/2 zz-rotation is then approximately 25 times 
the qubit period time (h/2Ej). 

Assuming that the qubits are operated at the degener- 
acy point, fluctuations in the biasing current will cause 
pure dcphasing. Nevertheless, the qubits will be decou- 
pled from current noise to first order at zero current bias. 
However, since relatively long periods of qubit coupling 
are needed to perform practical control operations, sup- 
pression of bias current fluctuations might be essential 
[5]. 

Finally it should be emphasized that although this 
paper has been concerned with the charge qubit limit 
Ej/Ec <C 1, the analysis and the design for bias-current- 
controlled qubit-qubit coupling is equally relevant in the 
region of Ej/Ec ~ 1, characterizing the " Quantronium" 
charge-phase qubit [3] . A higher Ej jEc ratio introduces 
more charge states and flattens the bands, making the 
system less sensitive to background charge fluctuations. 
The coupling of neighbouring qubits will however still be 
controllable, and higher levels will not be excited during 
two-qubit gate operations provided that the bias current 
is switched on adiabatically. 

In conclusion, the present scheme provides a realis- 
tic solution for easy local control of the physical cou- 
pling of charge qubits via current biasing of coupling 
Josephson junctions or, alternatively, pairs of readout 
junctions. The design is in line with experimental main- 
stream development of charge qubit circuits and can eas- 
ily be fabricated and tested experimentally. Most im- 
portantly, it allows readout via currently tested meth- 
ods that promise single-shot projective measurement 
and even non-destructive measurements, via e.g. RF- 
reflection readout of a JJ threshold detector [25,26] or 
an SET [7,27]. The tunable coupling of the qubit chain 
allows easy implementation of CNOT and CNOT-SWAP 
operations. Independent two-qubit operations can be 
performed in parallel when the network consists of five 
qubits or more, and generalization to single-shot N-qubit 
gates seems possible. This may offer interesting new op- 
portunities for operating qubit clusters in parallel and 
swapping and teleporting qubits along the chain, for ex- 
perimental implementations of elementary quantum in- 
formation processing [20,31-35]. 
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